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Abstract  In  Monterey  Bay,  CA,  during  northwesterly,  upwelling  favorable  winds,  the  development  of  a 
southward  flowing  cold  jet  along  the  entrance  to  the  Bay  is  often  observed.  This  dense  cold  jet  separates 
warm  waters  of  the  anticydonic  circulation  offshore  from  the  water  masses  inside  the  Bay.  Interactions 
between  the  cold  jet  and  the  offshore  anticydonic  circulation  generate  ageostrophic  secondary  circulation 
(ASC)  cells  due  to  submesoscale  processes  as,  for  example,  flow  interaction  with  the  development  of  surface 
frontogenesis  and  nonlinear  Ekman  pumping.  Based  on  observations  and  modeling  studies,  we  evaluate 
the  impact  of  these  submesoscale  processes  on  the  formation  of  chlorophyll  a  filaments  during  late  spring- 
earlier  summer,  and  late  summer  time  frames.  We  show  that  during  the  late  summertime  frame,  ASC  leads 
to  the  development  of  filaments  with  high  values  of  chlorophyll  a  concentration  along  the  edge  of  the  cold 
jet-in  contrast  to  the  earlier  summer  time,  when  the  ASC  mixes  phytoplankton  much  deeper  to  the  area 
below  of  the  euphotic  depth,  and  chlorophyll  a  filaments  are  3-4  times  weaker. 
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1.  Introduction 

Recently,  there  has  been  an  increase  in  observational  and  modeling  studies  focusing  on  the  impact  of  sub¬ 
mesoscale  processes  on  biological  dynamics  in  the  ocean  [see,  e.g.  Levy  et  al.,  201 2,  review].  It  has  been 
found  that  submesoscale  processes  with  scales  of  0(1-10  km)  in  the  horizontal,  0(100  m)  vertical,  and  0(1 
day)  temporal  domains  impact  phytoplankton  growth  and  increase  productivity  in  a  variety  of  ways,  as  for 
example,  by  modulating  the  vertical  supply  of  nutrients  into  the  euphotic,  lighted  layer,  or  by  changing  the 
light  exposure  of  phytoplankton  by  modulating  the  strength  of  vertical  mixing  [Levy  et  al.,  201 2;  McGilli- 
cuddy  et  al.,  2003;  Gruber  et  al.,  201 1 ;  Capet  et  al.,  2008;  Mahadevan  and  Archer,  2000;  Allen  et  al.,  2005; 
Lapeyre  and  Klein,  2006;  Nagai  et  al.,  2008;  Calil  and  Richards,  2010;  Pall'as-Sanz  et  al.,  2010].  These  processes 
contribute  to  the  development  of  5-10  km  wide  submesoscale  filaments  in  the  form  of  streaky  chlorophyll 
patches  often  observed  in  satellite  images  of  ocean  color.  As  stated  in  Levyetal.  [2012],  the  impact  of  sub¬ 
mesoscale  processes  on  the  intensification  of  phytoplankton  filaments  depends  on  the  local  hydrography, 
euphotic  depth,  and  nutrient  distributions. 

During  upwelling  favorable  winds,  pockets  of  cold,  salty  water  appear  to  the  north  of  Monterey  Bay,  CA. 
These  pockets  grow  into  a  cold,  dense  jet  extending  offshore  and  flowing  southward  along  the  entrance  to 
the  Bay  [Rosenfeld  etal.,  1994;  Ramp  eta!.,  2009].  The  cold,  dense  jet  interacts  with  the  warmer,  less  saline 
anticydonic  circulation  offshore  [Breaker  and  Broenkow,  1 994;  Rosenfeld  et  al.,  1 994;  Paduan  and  Rosenfeld, 

1 996;  Ramp  et  at.  2009].  This  feature  is  associated  with  the  larger-scale,  meandering  California  Current  Sys¬ 
tem  and  is  not  locally  generated  [Brink  et  al.,  1991;  Breaker  and  Broenkow,  1 994;  Rosenfeld  et  al.,  1 994].  The 
anticydonic  California  Current  meander,  also  some-times  referred  to  as  the  Monterey  Bay  Eddy  (MBE),  is  a 
frequently  observed  feature  of  the  region  during  the  upwelling  [Rosenfeld  etal.,  1994;  Ramp  etal.,  2005, 
2009].  The  MODIS-Aqua  SSTsfor  11,13,  and  15  August  2003  (Figure  la)  illustrate  well  the  described  dynam¬ 
ics  during  the  upwelling  event:  the  development  and  intensification  of  a  southward  flowing  cold  jet  and  its 
interaction  with  the  warm  offshore  anticydonic  circulation.  The  MODIS-Aqua  chlorophyll  a  images  show 
streaky  filaments  (5-10  km  wide)  along  the  offshore  edges  of  the  jet  (Figure  1).  Figure  1  shows  that  fila¬ 
ments  are  connected  to  the  high  chlorophyll  water  masses  to  the  north  of  the  Bay,  and  because  of  a  pre¬ 
dominately  southward  flow  during  the  upwelling,  we  can  suspect  that  the  advection  of  relatively  high 
chlorophyll  water  masses  from  the  north  contributed  to  the  development  of  the  narrow  filaments  in 
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Figure  1.  MODIS  Aqua  SST and  Chlorophyll  a  images  for  August  2003  (a)  and  June  2008  (b).  Figure  la  Black  lines  on  MODIS  SST  and  Chlorophyll  a  images  indicate  AUV  DORADO  sections 
taken  11  and  1 6  August.  Figure  1b  Black  dots  on  MODIS  Chlorophyll  a  image  for  10  June  2008  indicate  locations  of  The  R/V  Point  Sur  stations. 


Figure  1.  At  the  same  time.  Figure  1  shows  that  these  narrow  filaments  are  maintained  in  the  frontal  area 
between  the  meandering  warm  anticyclonic  circulation  offshore  and  the  cold  jet  It  is  known[Co///  and  Rich¬ 
ards,  2010  ■,  Pallas-Sanz  et  al„  2010;  Levy  eta!.,  2012],  that  the  area  between  a  cold  jet  and  warm  anticyclonic 
circulation  is  prone  to  generation  of  submesoscale  processes,  which  contribute  to  the  development  of  nar¬ 
row,  streaky  chlorophyll  filaments  like  those  presented  in  Figure  1. 

It  is  challenging  to  reliably  quantify  the  contributions  of  phytoplankton  advection  versus  in  situ  proc¬ 
esses  to  the  development  and  intensification  of  submesoscale  filaments.  There  is  a  lack  of  high- 
resolution  surface  and  subsurface  (Lagrangian  and  Eulerian-type)  bio-optical  and  physical  observations 
and  the  present  maturity  level  of  bio-optical  models  does  not  provide  the  needed  capabilities,  especially 
at  submesoscale. 

In  the  present  paper,  using  available  bio-optical,  physical  observations  and  a  well-established  physical 
model,  we  want  to  address  the  following  questions:  Are  the  observed  streaky  chlorophyll  a  offshore  fila¬ 
ments  (Figure  1 )  in  the  Monterey  Bay  area  locally  intensified  by  the  submesoscale  processes?  Observations 
of  ocean  color  show  different  intensities  of  observed  submesoscale  filaments  (for  example,  the  intensity  of 
the  August  2003  filament  is  3-4  times  higher  than  the  June  2008  filament,  Figure  1).  What  are  the  reasons 
for  the  observed  differences  in  chlorophyll  a  concentrations  of  the  August  2003  and  June  2008  filaments 
(Figure  1)?. 

To  answer  these  questions,  we  have  organized  the  paper  as  follows:  Section  2  describes  observations  and 
models  used  in  this  study.  Results  are  presented  in  section  3.  Section  4  is  devoted  to  discussions  and 
conclusions. 

2.  Methods 

2.1.  Observations 

2.1.1.  Satellite  MODIS-Aqua  Ocean  Color  Data 

The  MODIS-Aqua  satellite  imagery  was  processed  using  the  NRL  Automated  Processing  System  (APS).  APS 
is  a  complete  end-to-end  system  that  includes  sensor  calibration,  atmospheric  correction  (with  near- 
infrared  correction  for  coastal  waters),  and  bio-optical  inversion.  APS  incorporates,  and  is  consistent  with, 
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the  latest  NASA  MODIS  code  (SeaDAS)  [Gould  et  a!.,  201 1 ;  Martinolich  and  Scardino,  201 1  ].  In  this  study,  esti¬ 
mates  of  the  chlorophyll  a  (Chi)  from  MODIS-Aqua  imagery  for  11-15  August  of  2003  and  10  June  of  2008 
were  used.  Chlorophyll  data  are  derived  by  OC3M  algorithm  [ O'Reilly  et  at.,  2000]  at  1  km  pixel  resolution. 
The  diffuse  attenuation  coefficient  Kd  (488)  was  estimated  in  accord  with  Lee  et  at.  [2005].  The  Kd  (488)  was 
used  to  estimate  the  euphotic  depth  (noted  Eu);  the  Eu  was  estimated  as  the  depth  where  photosynthetic 
available  radiation  (PAR)  is  1%  of  its  surface  value  [ Lee  et  al.,  2007]: 


Eu=  In  (0.01  )/Kc/(488) 


2.1.2.  Mooring  Data 

Observations  of  winds  and  PAR  from  the  Monterey  Bay  Aquarium  Research  Institute  (MBARI)  surface  moor¬ 
ings  Ml  (122.02°  W,  36.74°  N)  and  M2  (1 22.40°  W,  36.67°  N)  are  used  in  this  study  (Figure  2).  Near-surface 
3m  wind  speed  and  direction  were  measured  by  a  MetSys  monitor.  The  observed  PAR  was  measured  by  the 
spectroradiometer  mounted  on  moorings  approximately  3  m  above  the  water  surface  [Chavez  et  al.,  2000], 

2.1 .3.  HF  Radar  Surface  Currents 

Surface  current  observations  used  in  this  study  were  derived  from  the  California  Coastal  Ocean  Current 
Mapping  Program's  HF  radar  network  (www.cocmp.org).  Surface  currents  were  estimated  based  on  inputs 
from  five  HF  radar  sites  for  August  2003  and  seven  HF  radar  sites  for  June  2008.  Vector  currents  were  esti¬ 
mated  on  a  Cartesian  grid  with  a  horizontal  resolution  of  3  km  by  computing  the  best-fit  vector  velocity 
components  using  all  radial  velocity  observations  within  a  radius  of  3  km  for  each  grid  point  each  hour 
[Paduan  et  al.,  2006]. 

2.1.4.  Water  Samples 

The  R/V  Point  Sur  occupied  hydrographic  and  optical  stations  from  2  to  1 3  June  2008.  Temperature  and 
salinity  depth  profiles  binned  to  1  m  vertical  resolution  were  derived  from  Sea-Bird  SBE  9+  CTD  meas¬ 
urements  using  standard  Sea-Bird  processing  software.  Nutrients  analyses  were  performed  on  water 
samples  from  the  Niskin  bottle  rosette  in  accord  with  Pennington  and  Chavez  [2000].  Water  samples 
(540  mL)  collected  from  near-surface  (~0.5m)  Niskin  bottles  were  filtered  onto  Whatman  glass  fiber  fil¬ 
ters  (GF/F).  High  performance  liquid  chromatography  (HPLC)  analysis  provided  pigment  indices  and  pro¬ 
portion  factor  for  microplankton,  nanoplankton,  and  picoplankton  [Vidussi  et  al.,  2001],  The  pigment 
data  indicated  that  the  microplankton  fraction  was  predominantly  composed  of  diatoms  (based  on  the 
presence  of  fucoxanthin). 

2.1.5.  Autonomous  Underwater  Vehicles  (AUVs) 

Propeller-driven  AUV  such  as  the  MBAFU  manufactured  DORADO  have  been  described  in  Bellingham  etal. 
[2000]  and  Ryan  etal.  [2009].  Deployments  of  the  DORADO  during  1 1  and  16  August  2003  are  used  in  this 
study.  The  instruments  on  board  included  CTD,  optical  backscattering,  chlorophyll  fluorescence,  nitrate,  and 
bioluminescence  potential  sensors. 

2.2.  Model 

The  Monterey  Bay  circulation  model  used  in  this  study  has  been  described  in  detail  elsewhere  [Shul- 
man  et  al.,  2010,  201 1].  The  model  domain  is  shown  on  Figure  2a.  The  model  is  based  on  the  Navy 
Coastal  Ocean  Model  (NCOM),  which  is  a  primitive-equation,  3-D,  hydrostatic  model.  It  uses  the 
Mellor-Yamada  level  2.5  turbulence  closure  scheme,  and  the  Smagorinsky  formulation  for  horizontal 
mixing.  The  Monterey  Bay  model  uses  the  Navy  Coupled  Ocean  Data  Assimilation  (NCODA)  system 
[Shulman  et  al.,  2011].  The  NCODA  assimilates  satellite  altimeter  observations,  satellite  sea  surface 
temperature,  as  well  as  available  in  situ  vertical  temperature  and  salinity  profiles  from  ships  and 
gliders.  The  Monterey  Bay  model  predictions  during  August  of  2003  and  June  2008  were  validated 
against  assimilated  observations  (to  verify  that  assimilation  decreases  misfits  between  the  model  and 
assimilated  observations),  as  well  as  against  nonassimilated  observations:  temperature,  salinity,  and 
currents  profiles  at  moorings,  HF  radar-derived  currents,  and  aircraft  SSTs  [Shulman  et  al.,  2009,  2010, 
2011,  2013].  The  Monterey  Bay  model  is  forced  with  surface  fluxes  from  the  Coupled  Ocean  and 
Atmospheric  Mesoscale  Prediction  System  (COAMPS)  [Doyle  et  al.,  2009]  at  3  km  horizontal  resolution. 
Comparisons  of  COAMPS  predictions  with  observations  from  moorings  and  aircraft  surveys  were 
reported  in  Ramp  et  al.  [2009],  Doyle  et  al.  [2009],  Chao  et  al.  [2009],  and  Shulman  et  al.  [2009, 

2010].  Comparisons  demonstrated  high  complex  correlations  between  the  observed  and  COAMPS- 
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Figure  2.  (a)  Model  domain  with  locations  of  moorings  Ml  and  M2,  (b)  Observed  wind  velocities  at  moorings  during  1 1  16  August  2003  and  (c)  5  10  June  2008.  (d)  Observed  HF  Radar 
surface  currents  averaged  overlS  17  August  2003  (d),  and  for  5  June  2008  (e). 

predicted  wind  velocities,  and  good  agreement  in  magnitudes  and  the  extent  of  each  observed 
upwelling  event.  Comparisons  of  observed  and  COAMPS-predicted  short-wave  radiations  (SWR)  dem¬ 
onstrated  the  overestimation  of  the  SWR  in  the  COAMPS  predictions  at  mooring  locations,  and  for 
this  reason  COAMPS  SWR  values  were  adjusted  to  match  observed  mean  daily  values  of  SWR. 

In  this  study,  a  depth  of  mixed  layer  (MLD)  is  derived  from  the  Monterey  Bay  model  and  observations.  In 
both  cases,  the  MLD  is  computed  as  the  depth  at  which  temperature  deviates  by  0.3°C  from  the  near¬ 
surface  (at  2  m  depth)  temperature. 
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3.  Results 


The  development  and  intensification  of  a  southward  flowing  cold  jet  and  its  interaction  with  the  warm  off¬ 
shore  anticyclonic  circulation  are  well  illustrated  by  MODIS-Aqua  SSTs  for  August  2003  on  Figure  la  and  by 
HF  Radar  surface  currents  averaged  over  15-1 7  August  2003  on  Figure  2d.  At  the  edges  between  the  cold 
jet  and  warm  anticyclonic  circulation,  there  is  a  development  of  a  narrow,  submesoscale  (~5-10  km  wide) 
filament  with  high  chlorophyll  a  concentration  (Figure  la). 

Because  the  location  of  the  submesoscale  filament  (Figure  la)  coincides  with  the  area  between  the  cold  jet 
and  the  anticyclonic  circulation  offshore,  the  submesoscale  process  due  to  their  flow  interaction  with  the 
development  of  the  surface  frontogenesis  is  considered  here  [ Hoskin ,  1 982;  Levy  et  aL,  201 2;  Calil  and  Rich¬ 
ards,  2010].  Because  the  upwelling  favorable  winds  have  been  blowing  along  the  jet  and  over  the  anticy¬ 
clonic  circulation  for  more  than  a  week  (Figure  2b),  the  second  submesoscale  process  considered  here  is 
due  to  the  forcing  interaction — nonlinear  Ekman  transport  [Thomas  and  Lee,  2005;  Levy  et  al„  201 2;  Calil  and 
Richards,  2010],  Both  submesoscale  processes  are  likely  to  contribute  to  the  development  and  intensifica¬ 
tion  of  the  submesoscale  filament  [Calil  and  Richards,  2010]. 


In  the  first  process,  the  interaction  between  the  dense  cold  jet  and  lighter,  warmer,  anticyclonic  circulation 
leads  to  an  ageostrophic  secondary  circulation  (ASC).  These  ASC  cells  are  generated  in  a  plane  perpendicu¬ 
lar  to  the  density  front  [Hoskins,  1 982],  which  are  upward  (upwelling)  on  the  light  (anticyclonic)  side  of  the 
front  and  downward  on  the  dense  side  of  the  front.  The  ASC  cells  lead  to  a  restratification  flow  from  the 
light  site  to  the  dense.  In  Hoskin  [1982],  the  vector  Qi  is  used  to  qualitatively  diagnose  ageostrophic  vertical 
motion  due  to  frontogenesis; 


/  dug  db  dvg  db  dug  db  dvg  db  \ 
V  dx  dx  dx  dy’  dy  dx  dy  dy ) 


where  ug  and  vg  are  the  horizontal  geostrophic  velocities,  and  ||  and  are  the  horizontal  gradients  of 
buoyancy: 

b=  gf—  (2) 

Po 


with  p‘  being  the  deviation  from  the  reference  density  p0,  and  g  is  the  acceleration  due  to  gravity.  The  local 
maximum  of  Qi  creates  a  vertical  circulation  with  associated  upwelling/downwelling  ACS  cells.  Vertical 
velocity  due  to  frontogenesis  (noted  here  by  Waf)  is  estimated  from  the  vorticity  equation  [Hoskin,  1982] 
with  boundary  condition  at  surface  Wat  0. 

The  second  process  is  the  intensification  of  ASC  cells  due  to  nonlinear  Ekman  transport  [Thomas  and  Lee, 
2005;  Levyet  aL,  201 2].  This  is  a  result  of  downfront  winds  blowing  in  the  direction  of  the  frontal  jet.  When  a 
wind  blows  down  a  front,  cross-front  advection  of  density  by  Ekman  flow  results  in  a  destabilization  of  the 
water  column,  which  results  in  a  convection  that  is  localized  to  the  front.  The  vertical  circulation  associated 
with  the  ASCs  is  characterized  by  subduction  on  the  dense  side  of  the  front  and  upwelling  along  the  frontal 
interface.  It  creates  upwelling  on  the  anticyclonic  circulation  site  (with  negative  relative  vorticity)  and 
downwelling  on  the  jet  site.  According  to  Thomas  and  Lee  [2005],  the  dominant  term  in  nonlinear  Ekman 
transport  is: 


Me= 


Ta 

P0(f+Z) 


(3) 


where  Me  is  nonlinear  Ekman  transport,  x„  is  along-front  wind  stress  component,  f  is  Coriolis  parameter, 
and  £  is  relative  vorticity  of  geostrophic  flow,  and  the  leading  order  term  of  the  nonlinear  Ekman  pumping 
(NLEP)  is: 


Wnkp=— z-  VX— (4) 
Po  (f+£) 

where,  x  is  surface  wind  stress. 

Our  objective  is  to  demonstrate  whether  surface  frontogenesis  (represented  by  QI  (1))  and  nonlinear 
Ekman  transport  (represented  by  Me  (3))  contributed  to  the  development  of  chlorophyll  a  filament 
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Figure  3.  (top)  the  model  SST  (15  August  2003)  with  the  location  of  crass  section;  (bottom)  the  model  vertical  velocity  with  the  model  potential  temperature  isotherms. 


(Figure  la).  The  local  maxima  of  Q,  and  Me  create  a  vertical  circulation  with  associated  upwelling/downwel- 
ling  ACS  cells.  For  both  processes,  ASC  cells  are  upward  (upwelling)  on  the  light  (anticydonic)  side  of  the 
front,  and  downward  on  the  dense  side  of  the  front.  The  ASC  cells  modulate  the  mixing  of  material  down  to 
the  mixed  layer  depth  (MLD),  provide  injections  and  cycling  of  nutrients,  as  well  as  the  cycling  of  phyto¬ 
plankton  in  the  subsurface  down  to  the  base  of  the  mixed  layer  depth.  All  these  directly  influence  the  inten¬ 
sification  and  development  of  submesoscale  phytoplankton  filaments. 

Figure  3  shows  the  model  vertical  velocity  and  the  potential  temperature  for  a  section  across  the  cold  fila¬ 
ment,  the  warm  offshore  anticydonic  circulation,  and  the  frontal  area  between  them. 

The  model  vertical  velocity  and  subsurface  temperature  show  the  development  of  the  ACS  cell,  with  the 
ascending  part  of  this  cell  coinciding  with  the  warm  part  of  the  front.  This  is  in  agreement  with  what  would 
be  expected  from  the  surface  frontogenesis  and  nonlinear  Ekman  transport. 

Figures  4a  and  5  show  observed  and  modeled-predicted  properties  plotted  along  the  section  crossing  the 
submesoscale  filament  On  Figure  4a,  the  vector  Q1  magnitude  (1),  the  value  of  nonlinear  Ekman  transport 
Me  (3)  and  the  NLEP  (Wnkp)  were  estimated  by  using  the  COAMPS  wind  stress  and  the  Monterey  Bay  model 
fields  of  temperature,  salinity,  and  velocity  averaged  over  a  period  of  48  h  centered  on  15  August  00Z  2003. 
Estimated  values  of  Q1,  Me,  and  Wn,ep  were  averaged  over  the  top  15  m  depth  (Figure  4a,  third  row),  which 
is  the  depth  of  model-predicted  MLD  (Figure  4a,  second  row).  On  Figure  5,  the  model  vertical  velocity  (W), 
the  Wniep,  and  the  vertical  velocity  due  to  frontogenesis  (Wat)  ate  plotted  along  the  same  section  as  on  Fig¬ 
ure  4a.  All  three  velocities  were  averaged  over  a  period  of  48  h  centered  on  15  August  OOZ  and  over  the  top 
15  m  depth.  Locations  of  the  local  maximum  ofQI  (1)  and  Me  (3)  coincide  with  the  location  of  the  filament. 
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Figure  4.  Observed  and  model  predicted  properties  plotted  along  sections  crossing  filaments  on  (a)  15  00Z  August  2003  and  (b)  10  00Z  June  2008.  Sections  are  indicated  by 
white  lines  on  MODIS  Aqua  SST  images  (top  row).  Plots  along  the  section  of:  (second  row)  MODIS  Agua  chlorophyll  a  (red)  and  SST  (light  blue),  euphotic  depth  Eu  (black)  derived 
from  the  MODIS  Aqua  observations,  and  the  mixed  layer  depth  MLD  (green)  derived  from  the  model;  (third  row)  MODIS  Agua  chlorophyll  a  (red),  magnitudes  of  vector  Q1  (blue), 
nonlinear  Ekman  transport  Me  (magenta),  and  nonlinear  Ekman  pumping  (black),  which  are  averaged  over  top  15m  depth;  (bottom  row)  MODIS  Agua  chlorophyll  a,  magni 
tudes  of  vector  Ql,  nonlinear  Ekman  transport  Me  and  nonlinear  Ekman  pumping  Wniep,  which  are  averaged  over  top  40  m  depth.  Note,  that  scales  of  Wniqt  for  2003  and  2008 
are  different. 
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Figure  5.  The  model  vertical  velocity  W  (blue),  the  nonlinear  Ekman  pumping  Wntep  (black),  and  the  vertical  velocity  due  to  frontogenesis 
Wat  (magenta)  are  plotted  along  the  same  section  as  on  Figure  4a.  All  three  velocities  are  averaged  over  a  period  of  48  h  centered  on  15 
August  00Z  and  over  the  top  15  m  depth. 

In  the  area  of  the  chlorophyll  a  filament,  all  three  vertical  velocities  show  upwelling  on  the  warm  side  of  the 
front  and  downwelling  on  the  cold  jet  side,  which  is  in  agreement  with  the  development  of  ACS  cells  due 
to  the  nonlinear  Ekman  transport  and  the  surface  frontogenesis.  The  same  order  of  W„iep  and  W ^  indicate 
that  both  submesoscale  processes  (the  nonlinear  Ekman  pumping  and  the  surface  frontogenesis)  are 
important  processes  in  the  upper  ocean  dynamics  during  the  formation  and  development  of  the  filament. 
With  the  MLD  around  15  m  and  vertical  velocities  of  order  10  4  10  3,  the  time  scales  of  considered  sub¬ 
mesoscale  processes  are  around  0.1 7-1 .7  days.  These  time  scales  are  in  the  ranges  of  observed  time  scales 
of  phytoplankton  growth  rates  [ Landry  eta!.,  2009;  Waterhouse  and  Welschmeyer,  1995],  With  spatial  scales 
of  10-50  km,  the  time  scale  for  horizontal  advection  by  typical  surface  velocities  of  order  02-0.4  m/s  (Figure 
2d)  is  around  0.3-3  days.  Therefore,  considered  submesoscale  processes,  horizontal  advection,  and  phyto¬ 
plankton  growth  rates  had  similar  time  scales  during  the  development  of  the  observed  chlorophyll  a  fila¬ 
ment.  Note,  that  the  vertically  averaged  model  vertical  velocity  (W)  is  about  an  order  smaller  than  vertically 
averaged  Wniep  and  Waf.  Calil  and  Richards  [2010]  also  reported  larger  values  of  vertically  averaged  W„iep  in 
comparison  to  the  vertically  averaged  model  vertical  velocity. 

Figure  6  presents  subsurface  properties  along  the  AUV  DORADO  surveys.  On  16  August,  the  figure  shows 
the  deepening  of  isotherms  offshore  at  the  end  of  the  section  where  the  section  reaches  the  warm,  anticy- 
donic  circulation  offshore  (see  position  of  the  section  over  15  August  MODIS  SSTs  on  Figure  1  of  the  paper, 
top).  Also,  there  is  a  flattening  of  isotherms  as  the  section  moves  from  the  warm  anticyclonic  circulation 
through  the  frontal  area.  This  is  consistent  with  the  surface  frontogenesis,  which  tends  to  flatten  the  isopyc- 
nals/isotherms  [see,  e.g.,  Klein  and  Lapeyre,  2009],  In  the  area  of  filament,  both  model-predicted  MLD 
(~15  m.  Figure  4a,  second  row)  and  observed  MLD  (~15  m.  Figure  6, 16  August)  are  shallower  then  the  esti¬ 
mated  euphotic  depth  (Eu)  in  and  around  the  filament  (Figure  4a,  second  row).  In  this  case,  the  developed 
ASC  cells  kept  phytoplankton  in  the  lighted  area  and  supported  photosynthesis.  Also,  the  ASC  cells  pro¬ 
vided  a  vertical  supply  of  nutrients  into  the  euphotic  area  and  facilitated  consumption  of  nutrients  (see 
uplifting  of  nutridine  in  the  area  of  the  filament  on  16  August,  Figure  6). 

In  the  area  of  the  filament,  the  observed  changes  in  subsurface  properties  show  the  shallowing  of  MLD, 
increase  in  consumption  of  nutrients,  and  increase  in  subsurface  values  of  chlorophyll,  optical  backscatter- 
ing  and  bioluminescence  potential  from  1 1  to  1 6  August  (Figure  6).  DORADO  surveys  (16  August)  show  a 
coincidence  of  high  values  in  observed  chlorophyll  a,  BL  potential,  and  optical  backscattering  in  the  area  of 
the  filament  Moline  etal.  [2009]  demonstrated  that  such  a  coincidence  of  high  values  in  considered  bio- 
optical  parameters  is  associated  with  the  presence  of  bioluminescent  species  of  dinoflagellates  in  the  north¬ 
ern  part  of  the  Monterey  Bay.  Based  on  this,  we  suspect  that  the  observed  coincidence  of  high  values  of 
chlorophyll  a,  BL  potential,  and  optical  backscattering  point  to  a  presence  of  the  bioluminescent  dinoflagel¬ 
lates  in  the  offshore  submesoscale  filament. 
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Figure  6.  AUV  DORADO  surveys  taken  1  land  16  August  2003.  Isotherms  are  shown  in  black,  and  estimated  mixed  layer  depths  (MLD)  are  shown  in  magenta  on  temperature  plots  (top). 
Vertical  dashed  lines  indicate  location  of  the  high  chlorophyll  a  filament  1 1  August  where  the  filament  appeared  later  on  MODIS  Aqua  image  on  15  August  16  August  where  the  fila 
ment  observed  in  MODIS  Aqua  image  on  15  August. 


Figure  1  provides  comparisons  of  the  MODIS-Aqua  chlorophyll  a  and  SST  images  for  August  2003  and  June 
2008.  As  for  August  2003,  images  for  1 0  June  of  2008  also  show  a  development  of  a  chlorophyll  a  submeso- 
scale  filament  in  the  area  between  the  cold  dense  upwelling  jet  and  the  warm  anticydonic  circulation.  The 
horizontal  cross-shore  scales  of  both  (August  2003  and  June  2008)  filaments  are  very  similar  (5-1 0km).  How¬ 
ever,  the  maximum  of  the  June  filament  chlorophyll  a  is  about  3-4  times  smaller  than  in  August  2003  (Figure 
1).  Figure  1  and  2  indicate  that  hydrographic  conditions  were  very  similar  during  both  time  periods  (August 
2003  and  June  2008):  (1)  steady  upwelling  favorable  winds  with  comparable  magnitudes;  (2)  development  of 
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cold,  southward  flowing  jet  along  the  entrance  to  the  Bay,  and  (3)  the  presence  of  warm  anticyclonic  circula¬ 
tion  offshore.  Both  filaments,  August  2003  and  June  2008  (Figure  1),  coincide  with  areas  between  the  cold, 
dense  jet  and  the  warmer,  anticydonic  circulation  offshore.  This  indicates  that  ageostrophic  secondary  circula¬ 
tion  cells  due  to  submesoscale  processes  (surface  frontogenesis  and  nonlinear  Ekman  transport)  also 
impacted  the  development  and  intensity  of  the  June  filament  The  question  is  why  the  chlorophyll  a  concen¬ 
tration  of  the  June  2008  filament  is  3-4  times  lower  than  one  in  August  2003  (Figure  1)?. 

One  obvious  reason  is  the  temperature  difference  between  August  2003  and  June  2008  along  the  two 
patches.  Figure  4b  provides  plots  of  observed  and  model-derived  properties  for  10  June  of  2008  (which  are 
also  plotted  for  1 5  August  2003  on  Figure  4a).  The  difference  between  observed  in  August  and  June 
MODIS-Aqua  SSTs  is  ~3°C  (Figures  4a  and  4b).  The  impact  of  temperature  on  the  filaments  growth  can  be 
estimated,  for  example,  by  using  the  relation  between  the  maximum  of  phytoplankton  growth  and  the  tem¬ 
perature  T  [Eppley,  1972]: 

/w(7-)=tt,*1.066r  (5) 

where  is  the  maximum  growth  rate.  According  to  (5),  3°C  difference  will  result  in  a  ~21%  increase  in 
the  maximum  phytoplankton  growth  rate.  This  alone  does  not  explain  the  observed  3-4  times  difference  in 
intensity  of  phytoplankton  filaments. 

Table  1  presents  observed  nutrients  (nitrate  and  silicate)  at  three  stations  taken  on  10  June  2008  along  the 
weak  chlorophyll  a  filament  (Figure  1b).  According  to  Table  1,  nitrate  availability  on  June  10  2008  in  the  sub¬ 
surface  (down  to  40  m  depth)  is  much  higher  along  the  10  June  of  2008  filamentthan  along  the  AUV  DOR¬ 
ADO  survey  on  16  August  2003  (Figure  6),  indicating  higher  consumption  of  nitrate  in  August  of  2003  than  in 
June  of  2008.  HPLC  samples  analysis  (section  2)  demonstrated  that  microplankton  (composed  predominantly 
of  diatoms)  was  about  80-85%  of  total  phytoplankton  population  (diatoms  plus  nanoplankton  and  picoplank- 
ton)  in  the  top  0.5  m  at  three  stations  on  10  June  of  2008.  Therefore,  the  filament  of  10  June  2008  was  mostly 
composed  of  diatoms  (at  least  close  to  surface),  which  require  silicate(Si(OH)4)  for  growth.  According  to  Table 
1,  the  ratio  of  Si(0H)4/N03  is  around  1:1  or  higher  for  all  stations  and  depths,  indicating  ample  silica  was  avail¬ 
able  in  the  water  column  for  diatom  growth  [Brzezinski, 1985],  Also,  ratios  of  Si(0H)4/N03  at  stations  (Table  1) 
taken  along  the  filament  are  comparable  or  even  higher  than  observed  ratios  at  stations  (not  shown  here) 
taken  in  the  more  productive  northern  part  of  the  Bay  (Figure  1b).  All  the  above  indicate  that  conditions  were 
not  nutrient  limited  for  the  filament  intensification  and  development  in  June  2008. 

Observed  daily  averaged  PAR  values  at  mooring  Ml  were  around  140  W/m2  during  11-16  August  2003,  and 
during  5-10  June  2008.  At  mooring  M2,  the  observed  daily  averaged  PAR  values  were  higher  (~150  W/m2) 
during  5-10  June  2008  than  during  11-15  August  2003  (~130W/m2).  Therefore,  the  available  surface  light 
for  photosynthesis  was  comparable  for  both  time  periods.  Estimates  of  euphotic  depths  (Eu)  are  ~ 20-30  m 
for  15  August  2003  (Figure  4a,  second  row)  and  ~30  m  for  10  June  2008  (Figure  4b,  second  row)  in  the 
areas  of  two  filaments.  Therefore,  based  on  the  observed  PAR  and  estimated  euphotic  depths,  conditions  of 
light  availability  in  surface  and  subsurface  were  comparable  for  both  of  the  considered  time  periods. 

The  mixed  layer  depth  (MLD)  for  10  June  2008  was  estimated  using  the  same  criteria  as  for  1 5  August  2003 
(Section  2).  Based  on  temperature  presented  in  Tablel,  the  MLD  is  around  40  m  for  station  1  and  between 
40  and  60  m  for  stations  2  and  3.  Therefore,  MLD  for  10  June  2008  is  about  2-3  times  deeper  than  MLD 


Table  1 .  Observed  Temperature,  Nitrate,  and  Silicate  at  Hie  R/V  Point  Sur  Stations  Taken  on  10  June  2008 

T  (°Q  NOj  (/i  mole/kg)  SiO.  Or  mole/kg) 

Stations  Stations  Stations 


Depth 

1 

2 

3 

1 

2 

3 

1 

2 

3 

0 

9.98 

9.6 

9.96 

213 

23.0 

21.4 

263 

229 

207 

5 

997 

959 

9.93 

23.9 

24.3 

215 

279 

245 

21.1 

10 

997 

9.58 

9.93 

235 

24.6 

21.6 

27  2 

245 

21.1 

20 

992 

9.58 

9.91 

23.7 

24.9 

213 

286 

2SO 

210 

30 

986 

95 

9.86 

235 

22.9 

21.9 

286 

22.1 

217 

40 

966 

9.3S 

9.69 

24.6 

27.2 

22.7 

288 

288 

23.1 

60 

941 

9.03 

9.64 

25.0 

29.1 

23.4 

328 

338 

26.4 
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observed  during  August  2003  (~15  m),  and  deeper  about  10-20m  than  the  estimated  euphotic  depth  Eu 
(Figure  4b).  Bear  in  mind  that  the  Eu  for  August  2003  is  deeper  than  observed  MLD.  Therefore,  the  discussed 
above  ACS  cells  in  June  of  2008  (due  to  surface  frontogenesis  and  nonlinear  Ekman  transport)  transported 
phytoplankton  out  of  the  euphotic  zone,  reduced  the  presence  of  phytoplankton  in  the  lighted  area,  and 
reduced  the  exposure  of  phytoplankton  to  light,  and,  therefore  limited  growth  of  the  filament  in  June  2008. 
The  Q1  magnitude  (1),  Me  (3),  and  Wnkp  (4)  were  estimated  by  the  same  procedure  as  for  1 5  00Z  August 
2003.  Values  of  the  Q1  magnitude,  Me,  and  Wn(ep  (averaged  over  40m  (the  observed  MLD  for  June  2008)) 
are  larger  for  10  June  2008  (Figure  4b,  bottom  row)  than  corresponding  values  for  1 5  August  2003  (Figure 
4a,  bottom  row).  Note,  that  the  MLD  derived  from  the  model  predictions  for  10  June  of  2008  (Figure  4b,  sec¬ 
ond  row)  is  somewhat  shallower  than  the  observed  MLD  during  this  period  (based  on  observed  tempera¬ 
ture,  Table  1);  therefore,  we  should  expect  that  actual  values  of  Ql,  Me,  and  W^p  are  even  larger  than  the 
model  predicted  for  10  June  2008,  and  certainly  larger  than  the  values  for  1 5  August  2003.  Therefore,  ACS 
cells  during  June  2008  mixed  the  phytoplankton  and  nutrients  below  the  euphotic  depth,  while  providing  a 
slower  supply  of  nutrients  into  the  euphotic  depth  in  comparison  to  the  15  August  2003  conditions,  when 
the  MLD  was  shallower  than  the  euphotic  depth. 

4.  Discussions  and  Conclusions 

MODIS-Aqua  chlorophyll  a  imagery  has  demonstrated  the  development  of  submesoscale,  high  chlorophyll 
a  filaments  during  the  upwelling  events  of  August  2003  and  June  2008.  These  filaments  are  connected  to 
the  high  chlorophyll  water  masses  to  the  north  of  the  Bay  (Figure  1 ),  suggesting  that  the  advection  of  high 
phytoplankton  water  masses  from  the  north  might  have  contributed  to  the  development  of  the  observed 
submesoscale  filaments.  At  the  same  time,  these  filaments  are  confined  to  the  narrow  frontal  area  between 
the  cold  jet  and  the  meandering  warm  anticydonic  circulation  offshore.  This  frontal  area  is  prone  to  the 
development  of  the  submesoscale  processes  and  ASC  cells,  which  are  known  to  contribute  to  the  develop¬ 
ment  of  narrow,  submesoscale  chlorophyll  filaments  in  other  areas  of  the  world  [see,  e.g.,  Calil  and  Richards, 
2010;  Pallas-Sanz  et  a!.,  2010;  Levy  ef  a!.,  201 2],  As  we  have  stated  in  the  Introduction,  it  is  difficult  to  sepa¬ 
rate  the  effects  of  advection  versus  in  situ  processes  based  on  existing  observational  assets  and  bio-optical 
models,  especially  at  submesoscale.  For  this  reason,  the  focus  of  the  paper  is  on  whether  submesoscale 
local  processes  have  contributed  to  the  intensification  and  dynamics  of  these  observed  filaments. 

We  used  the  vector  Q\  and  the  vertical  velocity  due  to  frontogenesis  (Waf)  to  diagnose  ageostrophic  vertical 
motion  due  to  frontogenesis.  The  leading  order  terms  of  nonlinear  Ekman  transport  Me  and  the  nonlinear 
Ekman  pumping  W„/ep  were  used  to  qualitatively  diagnose  the  ageostrophic  vertical  motion  due  to  nonlin¬ 
ear  Ekman  transport. 

Plots  of  the  observed  and  model-predicted  properties  across-filament  (1 5  August  2003)  have  shown  loca¬ 
tions  of  the  local  maximum  of  Ql  (1)  and  the  nonlinear  Ekman  transport  Me  (3)  in  the  area  of  the  observed 
submesoscale  filament  (Figure  4a).  The  nonlinear  Ekman  pumping  and  the  vertical  velocity  due  to  fronto¬ 
genesis  show  upwelling  on  the  warm  side  of  the  front  and  downwelling  on  the  cold  jet  side,  which  is  in 
agreement  with  the  development  of  ACS  cells  due  to  the  nonlinear  Ekman  transport  and  the  surface  fronto¬ 
genesis.  The  same  order  of  W„/ep  and  Waf  indicate  that  both  submesoscale  processes  (the  nonlinear  Ekman 
pumping  and  the  surface  frontogenesis)  were  equally  important  during  the  formation  and  development  of 
the  filament  During  the  development  of  the  observed  chlorophyll  a  filament  both  submesoscale  processes 
had  comparable  time  scales  with  the  time  scales  of  horizontal  advection  and  with  time  scales  of  phyto¬ 
plankton  growth  rates  documented  in  the  literature. 

The  observed  temperature  from  the  AUV  DORADO  survey  has  shown  a  flattening  of  isotherms  as  the  sec¬ 
tion  moves  from  the  warm  anticydonic  circulation  through  the  frontal  area.  This  is  consistent  with  the  sur¬ 
face  frontogenesis,  which  tends  to  flatten  the  isopycnals/isotherms  [see,  e.g.,  Klein  andLapeyre,  2009], 

Because  the  observed  and  model-predicted  MLD  are  shallower  than  the  estimated  euphotic  depth  in  and 
around  the  filament,  the  developed  ASC  cells  kept  phytoplankton  in  the  lighted  area  and  supported  photo¬ 
synthesis.  Also,  the  ASC  cells  provided  a  vertical  supply  of  nutrients  into  the  euphotic  area  and  facilitated 
consumption  of  nutrients  (what  was  also  observed  by  the  DORADO  surveys).  All  the  above  indicate  that 
local  ASC  processes  maintained  phytoplankton  population  in  the  lighted  area,  supplied  nutrient-rich  subsur¬ 
face  waters  into  the  euphotic  area,  and  therefore  facilitated  the  intensification  and  development  of  the 
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observed  phytoplankton  rich  filament  during  1 1-16  August  2003  (Figure  1 ).  This  provides  the  answer  to  the 
first  question  stated  in  the  Introduction  section  of  the  paper. 

Comparisons  of  chlorophyll  a  submesoscale  filaments  on  15  August  2003  and  10  June  of  2008  have  shown 
that  they  have  similar  horizontal  across-filament  scales.  However,  the  maximum  of  the  June  filament  chloro¬ 
phyll  a  is  about  3-4  times  smaller  than  in  August  2003  (Figure  1).  The  difference  between  the  observed 
MODIS-Aqua  SSTs  in  August  and  June  does  not  alone  explain  the  observed  3-4  magnitude  difference  in  the 
intensity  of  phytoplankton  filaments.  Observations  have  shown  that  hydrographic  conditions  and  light 
availability  in  surface  and  subsurface  were  very  similar  during  both  of  the  considered  time  periods  (August 
2003  and  June  2008).  An  analysis  of  water  samples  collected  along  the  June  2008  filament  has  shown  that 
conditions  were  not  nutrient  (nitrate  and  silicate)  limited  for  the  phytoplankton  growth.  The  mixed  layer 
depth  (MLD)  along  the  June  filament  (~4Q-60  m)  was  about  2-3  times  deeper  than  the  MLD  during  August 
2003  (~1 5  m),  and  about  10-20  m  deeper  than  the  estimated  euphotic  depth.  Therefore,  the  discussed 
above  submesoscale  processes  and  associated  ACS  cells  mixed  the  phytoplankton  and  nutrients  below  the 
euphotic  depth  in  June  2008,  which  provided  a  slower  supply  of  nutrients  into  the  euphotic  depth  in  com¬ 
parison  to  the  August  2003  conditions,  when  the  MLD  was  shallower  than  the  euphotic  depth.  During  June 
2008,  ACS  cells  reduced  the  presence  of  phytoplankton  in  the  lighted  area,  and  therefore,  reduced  exposure 
of  phytoplankton  to  light,  and  limited  the  growth  of  the  filament  in  comparison  to  August  2003. 

Our  results  show  that  during  the  late  summer  time  frame,  ASC  leads  to  the  development  of  phytoplankton- 
rich  submesoscale  filaments  offshore,  along  the  edge  of  the  cold  jet  (when  the  MLD  and  euphotic  depth 
are  comparable) — in  contrast  to  the  earlier  summer  time,  when  the  ASC  mixes  phytoplankton  much  more 
deeply  to  the  depth  of  the  MLD,  which  is  below  the  euphotic  depth,  and  submesoscale  chlorophyll  a 
patches  are  3-4  times  weaker. 

Our  results  illustrate  that  during  persistent  upwelling  favorable  winds,  submesoscale  processes  and  ACS 
cells  can  modulate  the  development  and  intensification  of  offshore  narrow  (5-10  km  wide)  phytoplankton 
filaments.  These  submesoscale  and  ACS  cells  can  incubate  the  phytoplankton  population  offshore  (as  for 
example,  bioluminescent  dinoflagellates  during  August  2003).  These  offshore  phytoplankton  filaments  can 
migrate  onshore  during  relaxed  winds  following  the  upwelling,  and  provide  an  additional  source  of  phyto¬ 
plankton  bloom  development  in  and  around  Monterey  Bay.  Therefore,  the  discussed  offshore  phytoplank¬ 
ton  filaments  may  be  a  factor  in  the  Bay  ecosystem  health,  as  for  example,  in  the  development  of  such 
events  as  harmful  algae  blooms  (HABs).  All  these  emphasize  the  importance  of  further  observational  and 
modeling  studies  of  these  submesoscale  processes  which  impact  the  development  and  intensification  of 
offshore  submesoscale  phytoplankton  filaments. 
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